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Abstract 


Owing to the broad energy coverage of Insight-HXMT in the hard X-ray band, we detected the highest energy of 
pulsation exceeding 200 keV around the 2017-2018 outburst peak of the first Galactic pulsating ultraluminous 
X-ray source (PULX) Swift JO243.6+6124, which is the highest energy detected from PULXs to date. We also 
obtained the highest energy of pulsation of every exposure during the outburst in 2017-2018, and found the 
highest energy is roughly positively correlated with luminosity. Using our newly developed method, we identified 


the critical luminosity being 4 x 10°8 ergs 


' when the main peaks of the low and high energy pulse profiles 


became aligned, which separates the fan-beam dominated and pencil-beam dominated accretion regimes. Above 
the critical luminosity, the phase of the main peak shifted gradually from 0.5 to 0.8 until the outburst peak in all 
energy bands is reached, which is in agreement with the phase shift found previously at low energies. Our result is 


consistent with what is derived from spectral analysis. 
Key words: (stars:) pulsars: individual (Swift JO243.6+6124) 


1. Introduction 


BeXRBs are a sub-class of high mass X-ray binaries 
(HMXBs) consisting of a neutron star (NS) accreting material 
and angular momentum from a Be companion star. They are 
among the brightest X-ray sources in the sky. Giant (Type II) 
outbursts are usually observed from BeXRBs which can exceed 
an X-ray luminosity of 10°” ergs~' and last for several weeks 
to even months. The giant outbursts provide a good opportunity 
to study the accretion process and emission properties of 
BeXRBs (Reig 2011). 

The radiation structure and properties of BeXRBs depend on 
the luminosity and accretion rate. At lower luminosities, the 
radiation is produced near the polar caps of the NS and escape 
mostly along the magnetic axis forming a “pencil” beam 
pattern (Basko & Sunyaev 1975; Burnard et al. 1991; Nelson 
et al. 1993; Becker et al. 2012), while at higher luminosities, 
materials are decelerated by the radiation-dominated shock and 
accretion columns are formed at the magnetic poles. In this 
case, the X-ray emission runs away predominantly through the 
column wall forming a “fan” beam pattern (Davidson 1973; 
Basko & Sunyaev 1976; Becker et al. 2012; Mushtukov et al. 
2015). The transition luminosity between these two accretion 
regimes is called the “critical luminosity” (L,,i,), the detection 
of which based on the different beam patterns was proposed by 


— X-rays: binaries — stars: neutron 


Gnedin & Sunyaev (1973). In the radiation model of Becker 
et al. (2012), when the luminosity is below L,,;, (subcritical 
accretion), pencil beam radiation dominates (mixture of pencil 
and fan beam), and above Lerit (supercritical accretion) fan 
beam radiation dominates. The observed pulse profiles depend 
on luminosity, energy and the geometry of the rotating NS in 
the observer’s reference frame. Many studies have shown that 
reflection by the NS surface and gravitational lensing effects 
result in more complex beam patterns than the simplified 
“pencil-fan” picture and thus more complex pulse profiles (see 
e.g., Poutanen et al. 2013; Lutovinov et al. 2015; Postnov et al. 
2015; Mushtukov et al. 2018; Markozov & Mushtukov 2024), 
as usually observed. 

Determining Lerit is not only crucial to identify the emission 
pattern transition which helps to probe the accretion state of 
BeXRBs, but also useful to estimate the magnetic field of the 
NS (see e.g., Wang & Frank 1981; Becker et al. 2012; 
Mushtukov et al. 2015). However, due to the complexity of the 
pulse profiles especially at low energies, the determination of 
Lorit using only low energy data (e.g., below about 30 keV) is 
highly uncertain. Hou et al. (2022) (and more details in the 
Appendix of this work) demonstrated that the beam pattern, 
thus Lerit, can be identified by judging whether the phases of the 
high and low energy main peaks are aligned in phase, since 
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high energy photons (e.g., >50 keV) preferentially escape from 
the wall of the column, i.e., in the fan beam pattern. 

Ultraluminous X-ray (ULX) sources are ultra-bright sources 
with apparent X-ray luminosity greater than 1 x 10°° ergs | 
(the Eddington limit of a ~10M. black hole). Known as 
pulsating ULXs (PULXs), some of ULXs are found to have 
X-ray pulsations, establishing that the central engines are NSs 
(Bachetti et al. 2014; Fürst et al. 2016; Israel et al. 
2017a, 2017b; Carpano et al. 2018; Chandra et al. 2020; 
Rodriguez Castillo et al. 2020; Vasilopoulos et al. 2020; Weng 
et al. 2017). See King et al. (2023) for a recent review on 
ULXs. Swift JO243.6+6124 is the first Galactic PULX with a 
peak luminosity of >1 x 10° ergs”! during its 2017-2018 
giant outburst. It has a Be companion star, a pulse period of 
9.8 s, an orbital period of 27.7 days and a distance of 6.8 kpc 
(See e.g., Tsygankov et al. 2018; Wilson-Hodge et al. 2018; 
Zhang et al. 2019; Doroshenko et al. 2020; Kong et al. 
2020, 2022; Wang et al. 2020). 

Insight-HXMT has a wide energy band of 1-250 keV, and 
the largest effective area above 80 keV. It is the only satellite 
that can measure pulse energies exceeding 100 keV and its 
wide energy range provides the opportunity to determine L,,it 
using the method that is described above. In previous studies, 
>130 keV pulsations have been detected from the PULX RX 
J0209.6-7427 (Hou et al. 2022) using Insight-HXMT data, and 
the method has been verified in two accreting X-ray pulsars 
EXO 2030+375 (Fu et al. 2023) and RX J0440.9+4431 (Li 
et al. 2023), respectively. 

In this paper, we search for the highest energy of pulsation of 
Swift J0243.6+6124 using the broad-band observations by 
Insight-HXMT, and determine Leit using our method. We 
describe the observation and data analysis in Section 2 and 
present our results in Section 3. We finally discuss our findings 
in Section 4 and summarize in Section 5. 


2. Observations and Data Analysis 


Insight-HXMT carries three slat-collimated X-ray tele- 
scopes: the High Energy X-ray telescope (HE, Liu et al. 
2020), the Medium Energy X-ray telescope (ME, Cao et al. 
2020), and the Low Energy X-ray telescope (LE, Chen et al. 
2020). Its broad energy band (1-250 keV), large detection area 
(5100 cm? in 20-250keV for HE) and the small dead time 
make it a powerful satellite in X-ray studies of bright X-ray 
sources (Zhang et al. 2020). 

Swift JO243.6+6124 was observed continuously by Insight- 
HXMT from 2017 October 7 (MJD 58033) to 2018 February 
21 (MJD 58170) covering its entire 2017-2018 outburst, 
providing rich data for investigating the properties of this 
source. Insight-HXMT made 102 pointing observations of 
Swift JO243.6+6124 during this outburst. The Insight-HXMT 
Data Analysis Software (HXMTDAS) v2.04 with default filters 
is used to reduce the data. These include screening data with 
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Figure 1. The luminosity evolution of Swift JO243.6+6124 as observed by 
Insight-HXMT from MJD 58033 to 58170 (2017 October 7 to 2018 February 
21) in the 2-150 keV range. The blue stars mark the nine exposures at the 
outburst peak. 


elevation angle (ELV)> 10°, geometric cutoff rigidity 
(COR) > 8 GeV, offset for the point position <0°04, and time 
beyond 300s to the South Atlantic Anomaly. The photon 
arrival time is corrected by the solar system barycenter 
(R. A. = 40°918, decl. = 61°4341) and by the binary orbit 
with the orbital parameters from Fermi/GBM measurements.’ 

To obtain phase-coherent pulse profiles, we fit spin 
frequency over time using the cubic spline interpolation, where 
the frequencies are from Fermi/GBM measurements, and then 
convert photon time after corrections to a sequential pulse 
phase (1), where ġo is set to 0 at tọ = 59027.49906951 (MJD), 
which is the epoch of the first Fermi/GBM_ periodicity 
detection (see Sugizaki et al. 2020, for details of this method). 
In addition, the dead time effect has been corrected for the 
pulse profiles of HE and ME using the recorded dead time in 
view of the high luminosity of Swift JO243.6+6124 and the 
corresponding significant dead time effect. 


3. Results 
3.1. Search for the Highest Energy of Pulsation 


Figure 1 shows the luminosity evolution of Swift JO243.6 
+6124 from MJD 58033 to 58170 (2017 October 7 to 2018 
February 21) in the 2-150 keV range, and the blue stars mark the 
nine exposures to be combined at the outburst peak. The 
luminosity was calculated by multiplying the flux by 4nd”, where 
the distance d was taken as 6.8 kpc, and the flux was estimated 
by fitting the broad-band Insight-HXMT spectra using con- 
stant*tbabs« (bbodyrad+Gausstcutoffpl) model. In 


7 https: //gamma-ray.msfc.nasa.gov /gbm/science/pulsars /lightcurves / 
swiftj0243.html 
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Figure 2. Pulse profiles in different energy bands by combining the nine 
exposures from 2017 November 3 to November 11 at the outburst peak of 
Swift J0243.6+6124. The pulse profiles are normalized by the mean count rate. 


order to search for the highest energy of pulsation, we combined 
nine exposures at the outburst peak between MJD 58060 and 
58068 (2017 November 3 to 2017 November 11, ObsID: 
PO11457701704-P011457702301). The pulse profiles of these 
nine exposures remained stable, therefore, it is reasonable to 
combine them. 

When folding pulse profiles using the events after barycentric 
correction, GTI boundary correction was performed, i.e., only the 
events in integer periods were adopted to fold the pulse profile. 
The pulse profiles of the nine exposures were created in different 
energy bands respectively and combined using the cross- 
correlation technique to align phases. Figure 2 shows the 
combined pulse profiles in different energy bands. We can see 
that the profile has one main peak and one minor peak which are 
separated by ~0.5 in phase. As the energy increases, both peaks 
retain until above 200 keV. The significance of the pulsation is 
quantified using the cross-correlation method that developed by 
Hou et al. (2022). We cross correlated the 200-240 keV profile 
with the profile in the lower energy band of 100-150 keV, which 
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Figure 3. Cross-correlation of the pulse profiles between 200-240 keV and 
100-150 keV. The upper panel shows the cross-correlation with phase lags. 
The lower panel shows the cross-correlation distribution from 10° Monte Carlo 
simulations of two non-correlated profiles. 


has similar profile shape but higher signal-to-noise ratio. The 
specific method is as follows: First, we calculate the cross 
correlation of the 200-240 keV profile with the 100-150 keV 
profile. Then we use the Monte Carlo method to simulate the 
cross-correlation distribution between two non-correlated profiles: 
one is generated by performing Poisson sampling of a flat profile 
and the mean of the simulated profile is taken as that of 
the 200-240 keV profile; the second profile is generated by 
performing Poisson sampling for each bin of the 100-150 keV 
profile. The cross-correlation of these two sampled profiles is then 
calculated, and finally the simulation is repeated for 10° times to 
obtain the distribution of the cross correlation. The result is shown 
in Figure 3. 

In the upper panel of Figure 3, the blue line indicates the 
phase lag at which the maximum cross-correlation value 
locates. It can be seen that the phase lag of the two profiles is 
small (~0.0451), which means that the profile of the high 
energy band is in phase with that of the lower energy band. In 
the lower panel, the blue line indicates the maximum cross- 
correlation value calculated from the upper panel. Then we 
evaluate the significance of such a cross correlation by 
calculating the proportion of the part on the right of the blue 
line, i.e., the chance probability (p-value) that the result is 
obtained by chance assuming that the null hypothesis (the 
200-240 keV profile is flat) is true. The significance is 4.260 
assuming a normal distribution. So we are confident that the 
highest energy of the pulsation of Swift JO243.6+6124 exceeds 
200 keV, which is the highest energy detected so far from 
known PULXs, exceeding the 130keV detected from RX 
JO209.6-7427 (Hou et al. 2022). 


Research in Astronomy and Astrophysics, 24:055006 (9pp), 2024 May 


° 
180 4 æ 
o 
e @e 
e eee o 
S 160 5 e e ooo 
g en G20 ot 
= e e o e e 
D o eo @ owo 
g 140 e coo o coco o oo 
€ eo o © @ecee o 
o e e omo 
S e œ 
= 1204 e o 
e e 
@e 
o 
100 4 e o oco 
e e oo e 
e e eo 
1038 1039 


Luminosity (2-150 keV, erg s~?) 


Figure 4. Evolution of the highest energy of pulsation with luminosity for all 
exposures from MJD 58033 to MJD 58170 covering the whole outburst of 
Swift JO243.6+6124. 


3.2. Evolution of the Highest Energy of Pulsation 


We also searched the highest energy of pulsation of every 
exposure of the entire outburst from 2017 October 7 (MJD 58033) 
to 2018 February 21 (MJD 58170) using the method described 
above, and investigated the evolution of the highest energy of 
pulsation with luminosity (2-150 keV). It appears that the highest 
energy of pulsation is roughly positively correlated with 
luminosity, as shown in Figure 4. To exclude the possibility that 
the evolution trend is entirely influenced by photon statistics, we 
randomly sampled 10° photons for different exposures to find 
the highest energy of pulsation. Figure 5 shows the sampled 
pulse profiles of two exposures at MJD 58043 and MJD 58065 
with luminosity of 3.7 x 10° erg s™! and 2x 10” erg a 
respectively. The highest energy of pulsation of MJD 58043 
reaches 110 keV, while that of MJD 58065 reaches 140 keV. 
Therefore, we propose that the highest energy of pulsation 
increases with the increase of the luminosity and such evolution is 
intrinsic rather than pure statistical effect. 


3.3. Determine the Critical Luminosity 


We found that once the accretion column is formed at 
luminosity around (from below to above) Leit high energy 
photons (e.g., >50 keV) preferentially escape from the wall of 
the column, i.e., in the fan beam pattern. More specifically, 
they escape from the bottom wall of the column, since the 
optical depth ratio (7/71, i.e., parallel/perpendicular to the 
magnetic field) increases rapidly with energy, and it is thus 
hard for photons to escape from both the top and the top wall of 
the accretion column (Hou et al. 2022, and more details in the 
Appendix of this work). This means that we can use the highest 
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energy pulse profile to pin-point the spin phase when the 
radiation from the wall of the accretion column is observed 
with the maximum flux, i.e., the phase of the fan beam 
maximum, even at luminosity below L,,i. The highest energy 
pulse profile alone cannot be used to identify if the accretion is 
below or above Leit since the radiation below about 30 keV 
dominates the total X-ray radiation and thus the radiation at 
higher energy is only a very small part of the total X-ray 
luminosity. In other words, we need to find out the transition 
luminosity when the low energy pulse is dominated by the fan 
beam emission. However, lower energy photons can contribute 
to both the pencil and fan beams, though predominately in the 
pencil beam below Lerit or in the fan beam above Lit. 
Therefore, our strategy is to use the highest energy pulse 
peak to first pin-point the spin phase when the radiation from 
the wall of the accretion column is observed with the maximum 
flux, i.e., to find the peak phase of the highest energy emission, 
and then examine at what luminosity the low energy pulse peak 
becomes aligned with the high energy pulse peak of the fan 
beam pattern. This luminosity should be Leri We thus 
generated the pulse profiles evolving with energy for every 
single exposure. In order to cover a wide range of luminosity, 
we selected 12 exposures during the outburst decay as 
indicated in Figure 6, and presented their pulse profiles 
evolving with energy in Figure 7. According to the pulse 
profile evolution from Panels (12) to (5), i.e., with the 
increasing luminosity, the main peak in the high energy band 
of 50-100 keV remains at around phase of 0.5 (red dotted line). 
In the pulse profiles of Panel (7) for MJD 58102 and after, the 
main peaks in high and low energy bands are not in phase, 
indicating that it is the subcritical regime with pencil beam 
emission dominating. Then in Panel (5) for MJD 58 099, the 
main peak in the lower energy band (1-4 keV) is also at the 
same phase of 0.5, i.e., the main peaks in high and lower 
energy bands are in phase. So, this epoch is in the supercritical 
regime with fan beam emission dominating. Thus, we propose 
that the transition from pencil-beam dominated to fan-beam 
dominated emission regimes occurred at the epoch of Panel (5), 
and Lert is calculated correspondingly as 4 x 10°% erg s™!. 


4. Discussion 


Several attempts have been carried out to determine Lerit 
basing on timing and spectral properties of Swift JO243.6 
+0124 during its 2017—2018 outburst. Its pulse profiles show 
complex evolution with significant luminosity and energy 
dependence. Using NICER and Fermi/GBM observations, 
Wilson-Hodge et al. (2018) showed that near a luminosity of 
1 x 10°% erg s~', the pulse profiles evolved from single-peaked 
to double-peaked, the pulsed fraction changed from decreasing 
to increasing with increasing luminosity, and the hardness 
ratios became softening above | x 10°8 erg s~". These beha- 
viors repeated at the same luminosity during the outburst rising 
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Figure 5. The pulse profiles of two exposures (MJD 58043 and MJD 58065) at luminosities of 3.7 x 10° erg s~' (left) and 2 x 10° erg s~' (right), respectively, 
randomly sampled with 10° photons. The highest energy of pulsation of MJD 58043 reaches 110 keV, while that of MJD 58065 reaches 140 keV. 
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Figure 6. The same light curve of Swift JO243.6+6124 as Figure 1. The blue squares mark the exposures the pulse profiles of which are shown in Figure 7, and the 
numbers on the right of the squares are serial numbers corresponding to those of the pulse profiles shown in Figure 7. 


and decaying parts, and this luminosity was interpreted as the 
Lait Which represents the transition of two accretion regimes 
between a Coulomb collisional stopping mechanism and a 
radiation-dominated stopping mechanism. At about the same 
luminosity of 1.5 x 10°% erg s~', Doroshenko et al. (2020) also 


found the change of pulse profile shape using data from 
Insight-HXMT, and also considered it as a transition between 
two accretion regimes. Besides, they found at higher luminosity 
of about 4.4 x 10° ergs ', there was a change in both pulse 
profile shape and power spectrum, which they associated with 
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Figure 7. Pulse profiles of Swift J0O243.6+6124 evolving with energy of 12 exposures during the outburst decay. The red dotted lines in Panel (5) mark the phase of 
the main peak in high energy. In other panels, lines are marked at the same phase as those in Panel (5). 
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the transition of the accretion disk from a gas pressure 
dominated state to radiation pressure dominated state. 

However, Kong et al. (2020) considered 4.4 x 10°8 ergs | as 
Lait instead at which the accretion regimes transferred based on 
the observed spectral property that the photon indices changed at 
this luminosity and above it had positive correlation with 
luminosity. Kong et al. (2022) obtained another value of Lerit as of 
2.9 x 10°8 erg s7! using the relation Leit = 1.5 x L37BIg/'5 
(Becker et al. 2012), where the NS magnetic field B was estimated 
based on the detection of the cyclotron resonance scattering 
feature (CRSF) around 120-146 keV from this source. Liu et al. 
(2022) found there was a 0.25 phase offset between the double- 
peak profiles at luminosities below 1.5 x 10°% ergs~' and above 
5 x 108 erg s_', and considered 1.5 x 10°8 erg s | as Lorit and the 
onset of the accretion column. The critical luminosity of 
4 x 10* ergs”! determined from our method in this work is 
close to that obtained by Kong et al. (2020) based on the evolution 
of photon indices with luminosity. 

From Figure 7, we can see at even higher luminosities 
(Panels (1)-(4)) than Lerit (Panel (5)), the main peak in the high 
energy band shifted from phase 0.5 to around 0.8. The phase 
shift was also noted in Doroshenko et al. (2020), but was 
obtained based on lower energy pulse profiles and they did not 
show the case for the pulsation at higher energies, which is 
indeed our new result. We note that at such high luminosities, 
the accretion flow covering the magnetosphere of an NS 
becomes optically thick. This can strongly influence the 
formation of pulse profiles resulting in, for example, increased 
pulse fraction and smooth profiles (Mushtukov et al. 
2017, 2019; Brice et al. 2023), the latter is indeed observed 
in our analysis (e.g., Panels (1)-(2)). 

This phase shift may be explained by the gravitational 
lensing effect which caused light bending for the photons 
from the fan-beam of the opposite accretion column. In 
general, with the growing of the accretion column, i.e., 
increasing luminosity, we will see more and more photons 
originating from the opposite accretion column and both the 
low and high energy profiles will be dominated by the fan- 
beam emission pattern from both the facing and opposite 
accretion columns. Detailed discussion of the pulse profile 
evolution far above Leri, by taking into account of the 
gravitational lensing effect and the accretion column 
geometry, is out of the scope of this work. 


5. Summary 


In summary, we analyzed the pulse profile of Swift JO243.6 
+0124 based on observations from Insight-HXMT, and 
obtained the highest energy of pulsation and the critical 
luminosity. The highest energy of pulsation reached above 
200 keV, the highest observed so far from pulsating ultra- 
luminous X-ray pulsars. The critical luminosity was determined 
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using our method which involves judging whether the main 
peaks of the pulse profiles in low and high energy bands are in 
phase. This method has also been applied successfully to two 
X-ray binary sources EXO 2030+375 and RX J0440.9+4431. 
The critical luminosity that we obtained was about 
4 x 10°% ergs™!, roughly consistent with that obtained from 
spectral analysis in previous works. 
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Appendix 


The critical luminosity distinguishes which one dominates 
between pencil-beam and fan-beam radiation patterns. The 
pencil beam is parallel to the magnetic field axis, while the fan 
beam is perpendicular to the magnetic field. The optical depth 
determines whether a photon can escape from the accretion 
column, and is defined as (Arons et al. 1987): 


T = ndl, (A1) 


where ne is the electron number density, o is the scattering 
cross section between photons and electrons, and / is the 
photon propagation length. It is theoretically possible to 
determine which radiation pattern dominates by comparing 
the optical depth parallel to and perpendicular to the magnetic 
field. Here we used the simplified Compton scattering cross 
sections, while the scattering cross sections in a strong field 
limit were calculated accurately by various authors (Her- 
old 1979; Daugherty & Harding 1986; Mushtukov et al. 2016). 
The simplified cross sections depend on the photon energy, the 
magnetic field, the polarization mode, and the photon direction 
momentum: 


LE<E, 
(a) O-Mode 
oj = or(E/E,)’, 
OL = OT, (A2) 
(b) X-Mode 
ox = or(E/E)?, (A3) 
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Figure A1. Assuming O-mode photons with redshift of 0, an NS with mass of 1.4M, and radius of 10 km, the ratio of 7) to 7, evolves with energy, magnetic field 


and luminosity. Only the cases with the ratio larger than 0.1 are shown. 


2. ESE, 
(A4) 


00 = Ox = OT, 


where the O-mode is the ordinary polarization mode, and its 
electric vector is in the plane composed of the magnetic field B 
and the photon momentum hk, the X-mode is the extraordinary 
polarization mode in which the electric vector is perpendicular 
to the plane of magnetic field and photon momentum, ør is the 
Thomson scattering cross section, and E, = 11.6B,2/(1 + 2) is 
the CRSF line energy, where Bj. is the magnetic field in units 
of 10'2 G. ne and l are calculated using the analytical 
expressions in Becker et al. (2012). We consider the case 
where there is radiation-dominated shock such that the 
accretion column is formed. n, can be derived from the mass 
conservation relation, M = TNeMpV, where v equals to the 
post-shock velocity of the infalling matter, which is 1/7 of 
the freefall velocity approaching the top of the shock, and ro is 


the radius of the accretion column. M is the accretion rate and 
is related to the X-ray luminosity as L = GMM /R. 

The photon propagation length / in the direction parallel to 
the magnetic field is estimated by the height of the radiation- 
dominated shock: 


-1 
H = 1.14 x 105 cm| —— (=) a | 
1.4 Mo 10 km /\ 10%” erg s7! 


(A5) 


and in the direction perpendicular to the magnetic field / is 
estimated by the accretion column radius rọ which is the 
maximum length in this direction: 


—1/2 —1/14 
r= 193 x 10%em( x 
0.1 1.4 Mo 


11/14 -2/7 1/7 
x( £ ) (at ) an E 
10 km 10!? G 10% erg s7! 
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Figure A2. Assuming O-mode photons with redshift of 0, an NS with mass of 
1.4Mo and radius of 10 km, the lines show the relation between the magnetic 
field and energy for different luminosities when the ratio of 7} to 7, is unity. 


In the calculation, we assume typical NS parameters: mass 
M = 1.4Mọ, radius R = 10 km, and redshift z= 0. 

For O-mode photons, we calculated the ratio of optical depth 
parallel and perpendicular to the magnetic field (7,/7,) for 
different photon energies, luminosities and magnetic field 
strengths, as shown in Appendix Figure Al. We can see that on 
the one hand, the ratio of optical depth increases rapidly with 
the increase of photon energy under different luminosities and 
magnetic field strengths, which suggests that the higher the 
energy, the more likely the photons escape from the wall of the 
accretion column in the form of a fan beam. On the other hand, 
the higher the luminosity, the larger the ratio, i.e., the easier the 
fan beam is dominated. 

In addition, Appendix Figure A2 shows the relation of the 
energies with the magnetic field strengths under different 
luminosities when the ratio equals unity. This defines the 
starting energy of the fan beam dominated pattern. We can 
draw two conclusions: (1) For a given luminosity, the stronger 
the magnetic field, the higher the starting energy of the fan 
beam dominated pattern; after crossing the starting energy, the 
higher the energy, the larger the ratio, the more that is 
dominated by the fan beam. (2) The higher the luminosity, the 
lower the starting energy of the fan beam dominated pattern, or 
the larger the energy range of the fanbeam dominated pattern, 
thus the easier it is to find the moment when the high and low 
energy pulses begin to be dominated by the fan beam pattern, 
i.e., they are in phase, and the luminosity at this moment will 
correspond to Lert. For X-mode photons, the scattering cross 
section is the same in all directions, such that the ratio of the 
optical depths only depends on the path ratio in both directions 
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and is independent of energy. Therefore, when the actual 
radiation is a mixture of photons of the two modes, the above 
conclusions are qualitatively unchanged. 
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